Optical Vortices with Starlight: Implications for ground-based stellar coronagraphy 
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Using an I = 1 blazed fork-hologram at the focal plane of the Asiago 122 cm telescope, we 
obtained optical vortices from the stellar system Rasalgethi (a Herculis) and from the single star 
Arcturus (a Bootis). We have analyzed the structure of the optical vortices obtained from non- 
monochromatic starlight under very poor seeing conditions using a fast CCD camera to obtain 
speckle patterns and carry out the lucky imaging technique, alternative to adaptive optics. With 
the insertion of a red filter and of a Lyot stop we performed 1 = 1 optical vortex coronography the 
double star HD74010. The results are in agreement with theory and numerical simulations. Our 
results open the way to applications of optical vortices to ground based astronomical observations, in 
particular for coronagraphy with I > 1 masks. No intrinsic orbital angular momentum was detected 
in the starlight. 
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Introduction - Optical vortices (OVs) are phase de- 
fects embedded in light beams endowed with orbital an- 
gular momentum (OAM). OVs, and more in general light 
beams carrying OAM, are generated after the insertion of 
a phase modulating device (PMD) that imprints a cer- 
tain vorticity on the phase distribution of the original 
beam. Such beams can be mathematically described by 
a superposition of Laguerre-Gaussian (L-G) modes char- 
acterized by the two integer-valued indices £ and p. The 
azimuthal index £ describes the number of twists of the 
helical wavefront and the radial index p gives the number 
of radial nodes of the mode. The electromagnetic field 
amplitude of a generic L-G mode in a plane orthogonal 
to the direction of propagation is 
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where wo is the beam radius and L™(x) is the associ- 
ated Laguerre polynomial. The presence of a phase fac- 
tor exp(— i£9) implies that these cylindrically symmetric 
modes carry an OAM equal to £h per photon, relative 
to their symmetry axis P, Q. For the same reason, a 
phase singularity is embedded in the wavefront, all along 
the propagation axis, with topological charge equal to £. 
The intensity distribution of an L-G mode with p = is 
generally shaped as a ring with a central dark hole, where 
the intensity is null due to total destructive interference. 
The radius of maximum intensity of this donut grows as 
the square root of £ and the intensity value decreases as 

Experimentally, these properties are produced with 
beams propagating through nonlinear optical systems Q 
and Kerr nonlinear refractive media [5]. This induced 
to argue that OAM could also be naturally generated by 
some astrophysical environments, possibly related to tur- 
bulent interstellar media with density discontinuities on 
wide scale ranges or to the distorted geometry around 



Kerr black holes. It was suggested that OAM could 
be also present in the blackbody radiation of the cos- 
mic microwave background . OVs have been already 
applied in diverse research fields such as laboratory op- 
tics, nanotechnologies and biology Q. Two astronomi- 
cal applications have been indicated: first, to overcome 
the instrumental limitations in the resolution of very 
close stellar sources due the Rayleigh separability cri- 
terion in diffraction-limited telescopes d, @, , second, 
to improve the capability of imaging extrasolar planets 
by peering into the darkness of an OV generated by a 
PMD inserted in the optical path of a Lyot-type corona- 
graph 0, [ll|, [l2| ■ Similar results can be obtained by di- 
rect imaging and analysis of tiny deviations of the main 
star's OV introduced by the presence of a close fainter 



companion [13j, [14[. However, ground-based telescopes 



will always feel the detrimental effects of atmospheric tur- 
bulence (seeing), which reduces the resolving power, even 
with adaptive optics. Therefore, it is important for the 
above astronomical applications of OVs to know how the 
seeing can affect the expected donut pattern generated by 
stellar sources. In this Letter we present the first images 
of stellar OVs generated by a PMD placed in the opti- 
cal path of the Asiago 122 cm telescope. In particular, 
we show how the effects of the atmospheric turbulence 
can be overcome to produce good quality OVs suitable 
for OV coronagraphy even under poor seeing conditions 
and in white light. Even if the contrast needed for the 
search of exoplanet or for classical coronagraphy can be 
achieved mainly with even £s , we decided to use I = 1 
OVs to avoid some experimental complications as higher 
order ones may not remain stable while propagating in 
turbulent media, and may split in first order ones j 1 61 ] . 

Non-monochromatic optical vortices - Several 
types of PMDs have been designed to generate axially 
symmetric OVs from an incident monochromatic on-axis 
beam. The most efficient are computer-generated fork- 
holograms [I3| and spiral phase plates (SPP) [HI]. How- 
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ever, the use of non-monochromatic beams in astronom- 
ical applications of OVs is mandatory to collect enough 
photons from faint stellar objects. In this case, the 
donut-sh&ped structure of monochromatic OVs will be 
changed pjJUj. Using SPPs, which are helicoidal trans- 
mission optical devices with a given total thickness vari- 
ation h s , all monochromatic OVs will have the same axis 
of symmetry. However, different wavelengths will present 
different OAM values according to £ = An(X) h s /X, 
where An is the difference between the refraction indices 
of the SPP material and the surrounding medium. Thus, 
the transmitted beam will possess different OAM val- 
ues producing <iorut£-shaped patterns of different sizes 
and a pro gressive filling of the central dark region. Re- 
cently 21| a scheme has been proposed to overcome this 
problem, but achromaticity is expected only for a limited 
bandwidth (~ 100 nm) in the visible. We have instead 
used a fork-hologram which is a grating with a number 
I of dislocations on its center. Monochromatic on-axis 
beams produce OVs with OAM indices £ — ml, where m 
is the diffraction order of the grating. Differently from 
SPPs, fork-holograms present the advantage of generat- 
ing OVs with the same £ at all wavelengths for on-axis 
polychromatic sources. For off-axis sources £ decreases 
as the star moves away from the center of the hologram. 

Spatial coherence of the incident beam is generally as- 
sumed to hold only for stellar sources when observed from 
space, but when ground-based telescopes are used, quite 
dramatic effects arc introduced by the atmospheric tur- 
bulence and spatial coherence is lost. However, for ex- 
posures shorter than the turbulence timescale (~ 10 — 
100 ms in the optical/near-infrared) we can obtain a 
group of bright speckles that represent the interference 
image produced by the coherent wavefronts generated by 
the random distribution of the atmospheric irregularities. 
A single nearly diffraction-limited stellar image can occa- 
sionally be produced when most of the stellar light falls 
in a single bright speckle. Fried [22I ] coined the term lucky 
exposures to describe high quality short exposures occur- 
ring in such a fortuitous way. This is the basis of lucky 
imaging [23| , one of the speckle imaging techniques used 
in modern astronomy as alternative to adaptive optics. 
Here we use the lucky imaging technique to obtain OVs 
from stellar sources useful for OV coronagraphy, with 
the additional restriction of selecting only those expo- 
sures where the star to be fainted is centered with the 
hologram dislocation. 

Telescope Observations - Here we report the results 
obtained in a campaign of observations carried out from 
2005 to 2007 with the Asiago 122 cm telescope, that fol- 
lowed some preliminary laboratory simulations [lif . We 
first observed the multiple system Rasalgethi (a Her) 
and the single star Arcturus (a Boo), a Her is a vi- 
sual binary composed by two unresolved binary systems 
presently separated by 4". 7: a Her A, formed by an 
M5 Ib-II semiregular variable (my = 2.7 — 4.0) and a 




FIG. 1: Optical setups, without (upper panel) and with (lower 
panel) spatial filter. T is the telescope; LI, L2, L3 the lenses; 
H the / = 1 fork-hologram; S the spatial filter slit. Stellar 
speckle patterns are sketched on the left of the optical setups, 
while the output images at the 0th and 1st diffraction order, 
on the right. All angles and displacements are exaggerated 
for clarity. 



fainter companion separated by 0".19 [24]], and a Her B, 
containing a GO II-III giant (my = 5.4) and a fainter 
secondary separated by 0".0035 25]. a Boo, instead, is 
a single star having visual magnitude my = 0.04 and 
spectral type K1.5 III. We also performed £ = 1 OV 
coronography with the double system HD74010, which 
F0 components are separated by 9". 9 and have visual 
magnitudes my^A — 7.72 and my t B = 7.87. The optical 
schemes for the two types of observations to image the 
OVs are sketched in Fig. [TJ the OV coronagraph is pre- 
sented in Our I = 1 fork- hologram H is blazed 
at the first diffraction order with 20 lines mm -1 and 
has an active area of 2.6 x 2.6 mm 2 . H was placed in 
proximity of the F/16 Cassegrain focus of the telescope. 
The OVs produced by non-monochromatic light beams 
crossing the fork-hologram show intensity patterns that 
appear as rings stretched along the dispersion direction 
with a central dark strip. The spectral dispersion also 
causes a partial filling of the central dark zone. Thus, 
if we want to use non-monochromatic OVs produced by 
fork-holograms for OV coronagraphy, we must limit the 
spectral range [13j and/or restore the donut shape (see 
e. g. [26}). In our case we used a variable spatial filter, 
S, made by a slit placed on the Fourier plane of the col- 
limating lens L2 to limit the dispersion of the light at 
the first diffraction order. This adjustable slit works as 
a tunable bandpass filter with flat spectral response and 
has been used only for single on-axis stars. We used a 
fast CCD camera with 660 x 494 and 7.4 x 7.4 fim pixels, 
16 bit dynamical range (4000 to 6700 A with a peak at 
5200 A) to image the OVs' shapes, and, for coronography, 
the ultrafast EMCCD Andor iXon+ camera. 

While observing the multiple system a Her we could 
simultaneously observe the zeroth and the first diffrac- 
tion order (£ — and £ = 1, respectively, for on-axis 
stars) . We can then verify how the observed OV patterns 
change with the distance of the beam from the center 
of the hologram. We placed the fork-hologram 30 mm 
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FIG. 2: Speckle pattern (right) and dispersed OVs (left) ob- 
tained by summing the 4 lucky exposures of a Herculis. The 
image is displayed in a squared greyscale. The reference sys- 
tem used for the extraction of the OV profiles is reported. 
The x axis corresponds to the direction of dispersion. 



before the telescope focal plane and set the main star 
(a Her A) at the center of the optical system. The mean 
seeing was ~ 3". 8. Under such bad seeing conditions the 
lucky imaging technique empirically predicts that the res- 
olution should improve by a factor 2 if we select 1% of 
frames taken with an exposure time of ~ 80 ms 23[ . We 
then recorded a sequence of 860 frames at a time step of 
70 ms. To select the best exposures we imposed the ad- 
ditional restriction that the intensity of the peaks of the 
a Her A OV were equal, that corresponded to the stel- 
lar beam intersecting the center of the hologram. Fig. [2] 
represents the sum of the 10% lucky speckle patterns ob- 
served together with their corresponding OVs. Since the 
signal from a Her A is saturated, an angular resolution 
of ~ 1".5 was derived from the full width at half maxi- 
mum (FWHM) of the intensity profile of a Her B. Being 
the central part of the OV profile of a Her A saturated, 
the choice of the lucky exposures has been made by ana- 
lyzing the intensity profiles taken across the wings of the 
dispersed OVs. In Fig. [3] we analyze in detail the OV 
profiles averaged over 10 pixels wide strips perpendicular 
to the dispersion, i. e. along the y axis of Fig. [5] We 
find that the ratio of the two maxima fitted with two 
Gaussians is Ra 1 = 0.995 ± 0.005 for the on-axis star 
and Rb> — 0.83 ± 0.01 for the off-axis star. We also nu- 
merically reconstructed the theoretical intensity profiles 
of the chromatically dispersed on-axis and off-axis OVs. 
These profiles present dips with intensity of 82% and 70% 
respectively. As expected for £ = 1, their central regions 
can never be totally obscured also under ideal seeing con- 
ditions. The profile A' of a Her A shows a central dip 
of 81%. The OV generated by a Her B, being off-axis, 
presents an asymmetric shape with a non-integer £ value. 

The single star a Boo was set at the center of the 
hologram and the spatial filter S was introduced to pro- 
duce a nearly-monochromatic circularly symmetric OV. 
We adopted a slit width of 0.1 mm, corresponding to a 
300 A bandpass width in the visible, that ensured enough 
S/N ratio for the 70 ms exposures. Since in this case 
the slit did not allow the simultaneous observation of 
the speckle patterns at the zeroth diffraction order, we 
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FIG. 3: Observed OV profiles of a Her A (A', left) and 
a Her B (£>' , right) obtained from the selected lucky expo- 
sures. The corresponding numerical simulations are repre- 
sented by the dotted lines. The thin lines are drown at any 
dip level. 



could only select a 17% sub-sample of frames that pre- 
sented OVs circularly symmetric and therefore centered 
with the optical singularity of the hologram. We assumed 
that this condition was met when the ratios of the two 
intensity peaks measured along two perpendicular axes 
across the OV were close to 1 within an error of ~ 10~ 2 . 
Fig. 0] shows the OV obtained by summing the selected 
lucky frames. The central region of the OV is not to- 
tally dark because of the loss of the starlight coherence 
due to the very bad seeing conditions and to the pres- 
ence of residual chromatism. We however notice that 
the contrast between the dark center and the bright ring 
improves down to 52% with respect to the previous un- 
filtered symmetric OVs. 

We finally observed the double system HD74010 
through a 100 A bandpass red filter centered at 6532 A 
without and with a Lyot stop. The exposure times were 
0.01 s. Fig. [5] shows the results of £ = 1 OV coronagra- 
phy: the intensity profiles of the two stars without and 
with the Lyot mask, obtained by averaging over a 40 pix- 
els strip, are shown in the left panels. Right panels show 
the corresponding snapshots. The on-axis component B 
of the binary system appears fainted by a factor ~ 1.7, 
close to the factor 2 derived from numerical simulations. 
The partial obscuration of the on-axis star is due to the 
fact that we were using a £ = 1 fork-hologram. Total 
obscuration may be mainly achieved with even £ values. 

Discussion and Conclusions - We demonstrated 
the feasibility of OV coronagraphy even under very bad 
seeing conditions with ground-based telescopes. Using 
the lucky imaging technique to improve the quality of our 
images, we obtained slightly symmetric OV patterns from 
on-axis starlight beams and performed a first trial of OV 
coronagraphy obtaining a fainting close to that that can 
be actually obtained with our I = 1 fork-hologram. We 
tested how the pattern of polychromatic starlight OVs 
depends from the off-axis position of the star beam. We 
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FIG. 4: Left panel: the OV obtained by summing the selected 
2% good frames (see text). Left panel: profile of the OV 
across the direction perpendicular to the dispersion. 
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FIG. 5: I = 1 OV coronagraphy of HD74010. Left panel: the 
averaged intensity profiles of the double star before (up) and 
after (down) the insertion of the Lyot stop. Right panel: the 
corresponding snapshots (see text). 



showed how polychromaticity and the lack of coherence 
produced by the atmospheric turbulence may heavily al- 
ter the theoretical donut profile of the OV reducing the 
contrast needed for coronagraphy. Although only even 
values of I can produce perfect rejection of the on-axis 
starlight as required for coronagraphy, we decided to start 
with an I = 1 fork-hologram to ensure the stability of the 
OVs along the optical path. This explains the impor- 
tance of this choice as a first test of OV coronagraphy. 
We found that our experimental results are consistent 
with theoretical predictions. We suggest that these new 
promising techniques could find their best application 
mainly at telescopes with adaptive optics, or in space 
instruments. Finally, since we did not detect any OV at 
the other diffraction orders of the hologram, we can con- 
clude that starlight does not have an appreciable intrinsic 
OAM. 

We would like to thank the Institut fur Expcrimcn- 
talphysik, University of Wien, Zeilingcr-Gruppe for sup- 
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